A detailed study of quasinormal frequencies of the Kerr black hole* 
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We compute the quasinormal frequencies of the Kerr black hole using a continued fraction method. 
The continued fraction method first proposed by Leaver is still the only known method stable and 
accurate for the numerical determination of the Kerr quasinormal frequencies. We numerically 
obtain not only the slowly but also the rapidly damped quasinormal frequencies and analyze the 
peculiar behavior of these frequencies at the Kerr limit. We also calculate the algebraically special 
frequency first identified by Chandrasekhar and confirm that it coincide with the n — 8 quasinormal 
frequency only at the Schwarzschild limit. 



I. INTRODUCTION 



The linear perturbations of the Schwarzschild spacetime were first examined by Regge and Wheeler 0, who were 
trying to solve the stability problem of the black hole. They derived what is now called the Regge- Wheeler equation, a 
. second order ordinary differential equation similar to the Schrodinger equation. The perturbation radiates away at the 
speed of light satisfying the wave equation and hence it corresponds to the gravitational wave around the black hole. 
Therefore to examine the perturbations of the black hole is equivalent to solving the scattering problem of gravitational 
t-H . waves around the black hole. Their pioneering work has led to the study of the quasinormal modes because the 
dominates most processes involving perturbed black holes. Quasinormal modes were first found by Vishveshwara 
and by Press g{ through numerical computations of the time evolution of the gravitational waves around the black 
hole. External perturbations excite the quasinormal modes which in turn appear as damped vibration, emitting 
gravitational waves with specific frequencies. This has been shown in several other investigations on the response 
■ of the black hole to the external perturbations. j| There have been several algorithms proposed to determine 
\Q \ these resonant frequencies of the black hole. The first method presented by Chandrasekhar and Detweiler J§| was not 
ON ■ sufficient to obtain accurate values of the quasinormal frequencies because their scheme, namely, a direct numerical 
"q ] integration of the Regge- Wheeler equation, was not stable, especially for rapidly damped modes. However, it is 
^1 remarkable that they were able to obtain for the first time values of the first few least damped quasinormal modes. 

Meanwhile, there have been some attempts to obtain the frequencies analytically. An inverse potential method by 
6J[)! Mashhoon g-JDJ was the first attempt to obtain the quasinormal frequencies using a potential for which the spectrum 
is analytically solvable. There is also a semianalytic approach, the WKB method used by Schutz and Will Jl2| , 
. \ and later refined by Iyer and coworkers JT||,[l4|]. They have provided relatively easy and intuitively understandable 
schemes for obtaining the quasinormal frequencies. The methods have also been applied to the Kerr and Reissner- 
Nordstrom black holes |lq] and have provided numerical values for fundamental modes. However, these algorithms 
cannot provide very accurate values for the frequencies and break down completely for rapidly damped modes. 

Leaver p7|| has shown that a continued fraction method previously used for determining the energy spectrum 
of the hydrogen molecule ion |l8| , |l{| can be generalized to determine the quasinormal frequencies. This method 
provides extremely accurate values for the frequencies of the Schwarzschild and Kerr black hole. It was found in 
the Schwarzschild case that there are an infinite number of damped modes for every multipole index I and these 
modes are located approximately parallel to the imaginary axis. Also for the Kerr black hole, he calculated how 
these modes move from the Schwarzschild quasinormal frequencies as the black hole parameter changes. Later, he 
generalized his method to the Reissner- Nordstrom black hole [|o|, and found out how the quasinormal modes move 
as the charge of the black hole changes. These results for the non-rotating cases are confirmed by other independent 
investigations: Nollert and Schmidt | pl| presented the Laplace transform method which are in excellent agreement 
with Leaver's results for the Schwarzschild black hole. Andersson and coworkers used the phase integral method 
for the Schwarzschild black holej22-p4] to obtain the consistent results. The last method was also used for the 
Reissner- Nordstrom black hole [E 5U26I1 . 
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There are several methods to precisely obtain the quasinormal frequencies for the non-rotating cases, nevertheless, 
the continued fraction method by Leaver is still the only known way which can be generalized to the Kerr quasinormal 
frequencies. The Kerr black hole is most important astrophysically, thus its quasinormal modes should be investigated 
in detail. At the same time, several questions were raised by the previous detailed studies of the quasinormal 
frequencies of the Reissner-Nordstrom black hole ]27j-p9t . In Ref. [E7| we have improved a continued fraction method 
for the extreme case and then discovered a curious fact that the gravitational and the electromagnetic frequencies are 
identical. Whether this is also the case for the Kerr black hole is still unknown. In Ref. we presented a detailed 
calculation of the quasinormal modes of charged black holes focusing on the rapidly damped modes for nearly extreme 
cases. We found that the rapidly damped modes show several peculiar features: the higher modes generally spiral 
into the value for the extreme black hole as the charge increases. We also found that there is a quasinormal mode 
that converges to the algebraically special mode |3(J when the charge vanishes, even though there was no correlation 
between these two modes in the charged case. It is still unresolved whether this is just a coincidence or can be 
explained physically. There are arguments for and against the existence of the quasinormal modes on the imaginary 
axis. Some methods indicate that there should be such a mode and that it corresponds to the ninth quasinormal 
mode. The problem is that the algebraically special mode of non-rotating black holes is purely imaginary and that 
most existing methods do not work for purely imaginary modes. In order to resolve this arguments, it will be useful 
to investigate the Kerr case because the algebraically special mode move away from the imaginary axis and thus 
can be compared with the quasinormal modes. It is also important to know in detail the behavior of quasinormal 
modes for nearly extreme Kerr black holes, because some astrophysically relevant black holes are a consequence of 
the coalescence of binary neutron stars and are believed to have a considerably large angular momentum. However, 
the behavior of these modes near the extreme limit is not yet well known. 

In this paper, the quasinormal frequencies of the Kerr black hole are numerically computed using Leaver's continued 
fraction method beyond his original calculation. Several unresolved problems about the Kerr quasinormal modes are 
discussed. 



II. CONTINUED FRACTION EQUATIONS 

It is well known that perturbations of the Kerr black hole are described by the Teukolsky equations. In case a t- 
and ip- dependence is given by e - lwt + JTOl P ; the separated differential equation for an angular part of perturbations is 

( —j— sixS^ 

W = a 2 w 2 u 2 - 2ausu + A lm + s - ■^ (2.2) 

1 — u z 



and that of a radial part is 



^-TT + 2 ( s + -M)^r + V(r)R s = 0, (2.3) 

K 2 — 2isK(r - M) 2 2 

V= —± + Aisujt - Ai m + 2aum- a z u z , (2.4) 



where 



u = cos 6>, (2.5) 
A = r 2 -2Mr + a 2 , (2.6) 
K = (r 2 + a 2 )u - am. (2.7) 

These equations correspond to the electromagnetic waves in the Kerr geometry when s = — 1 and the gravitational 
waves when s = —2. In what follows, we normalize 2M — 1 if otherwise mentioned. 



In the Schwarzschild case, Eq.(2.2) can be solved analytically and hence the separation constant Ai m can be obtained 
as Ai m — 1(1 + 1) — s (s + 1). However, in the Kerr case, we have to solve it numerically as follows Jl7| . Boundary 
conditions for Eq.(|2.2|) are that S s is regular at the regular singular points u = ±1. The indices there are determined 



by the regu larity at both singular points and are given by ±(m + s)/2 at U = 1 and ±(m — s)/2 at u = — 1. A solution 



to Eq.(2.2) can be expressed as 
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S s (u) = e- aau (l + u)l m - s l/ 2 (l - u )l™+ s l/ 2 On(l + u) n . (2.8) 

n=0 

The expansion coefficients are related by a three-term recurrence relation and the boundary condition at u = 1 is 
satisfied only by its minimal solution sequence. The recurrence relation is 

a ai + /3 a = 0, (2.9) 

a n a n+ i + (3 n a n + 7„a n _i =0, (n > 1), (2-10) 
where the recurrence coefficients are, with k± = |m — s\/2 and fc 2 = \m + s\/2, 

a„ = -2(n+l)(n + 2fci + l), (2.11) 
/3 n = n(n - 1) + 2n{k x + k 2 + 1 + 2aa) 

+2aa(2ki + s + 1) + (fci + /c 2 )(fci + fc 2 + 1) (2.12) 

-aV - s(s + 1) - ,4 im , (2.13) 

7„ = 2aLu(n + k 1 + k 2 + s). (2.14) 

The minimal solution sequence will be obtained if the angular separation constant A\ m is a root of the continue 
fraction equation below, 

„_a apji Q172 «273 «374 / 2 lg v 

h~ 02- 03 - 04-' 

The boundary conditions for quasinormal modes for the radial equation are 

R-2 ~ (r-r + )- s - la +e^ r , as r^r+, (2.16) 
R-2 ~ r - 1 - 2 * +<w e <wr ) as r^oo, (2.17) 



where + = (o;r+ — am)/b and & = vl — 4a 2 . Thus, the solution should be 

00 / 

R- 2 = e l " r (r-r-)- 1 - s+l " +la+ (r-r + y s - sa +^a n (— ?±\ , (2.18) 



n=0 x 

where the expansion coefficients are again defined by a three-term recurrence relation: 

a ai + /3 a = 0, (2-19) 

a n a n+ i + f3 n a n + 7„a„_i = 0, (n > 1). (2.20) 

The recursion coefficients are 

a„ = n 2 + (c + l)n + c , (2.21) 

/?„ = -2n 2 + (ci + 2)n + c 3 , (2.22) 

7„ = n 2 + (c 2 - 3)n + c 4 - c 2 + 2, (2.23) 



and c„ are defined by 



2i fu \ 

cq = 1 — s — iu) — — y — — amj , (2.24) 
Cl = -A + 2iiv(2 + b) + - am) , (2.25) 



b \2 
2i fui 



c 2 — s + 3 — Ziuj — — ^— — arri^ , (2.26) 
c 3 = Ll> 2 (4 + 2b- a 2 ) - 2amui - s - 1 + (2 + 



Auj + 2 i (uj 



b 



(f - ^) , (2.27) 



Aim + 

c 4 = s + 1 - 2cj 2 - (2s + 3)iw - — u [--am) . (2.28) 



4w + 2i ( u 

■ 1 * 1 



The boundary condition that Eq.(2.1S) converges as r — > 00 reduces again to 

n — B «07l QU72 «273 "374 

; ° 01- 02~ 03- 04-"' 

The problem is now to seek frequencies which satisfy both Eq.( 2.15| ) and Eq.( 2.29| ) 



«07l "172 "273 "374 ,„ 

= A) --5 — 75 — 75 — 75 (2-29) 

Pi— P2— P3~ P4 — 
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III. NUMERICAL RESULTS 



A. The case of m — 



The case of m = is similar to the Reissner-Nordstrom quasinormal modes and hence several things will be 
repeated. First, each to — mode always has a corresponding mode symmetric about the imaginary axis. This is due 
to a symmetry of Eqs. (2.15) and (2.29). It is obvious that these continued fraction equations are equivalent even 
after the transformation of to — ► —to, ilu — ► (iu>)*, and Ai m — > ^* TO - Therefore, when to = a solution always appear 
as a complex conjugate pair in (iu)). In contrast to the Reissner-Nordstrom black hole the mode- frequencies of the 
first few modes have the damping rates (imaginary parts of ui) which decrease monotonically as a/M increases. On 
the other hand, the oscillation frequencies (real parts of u>) increase with a/M, which is the same as the Reissner- 
Nordstrom black hole. These results are clear from Fig. 1, where we show the behavior of the first three modes for 
(— s, I, m) = (1,1,0), (2,2,0), (1,2,0), and (2,3,0). The left endpoint of each line corresponds to the Schwarzschild 
mode, and the right endpoint corresponds to the a/M = 0.99 Kerr mode. The electromagnetic modes of the Reissner- 
Nordstrom, which correspond to the s = — 1 modes here, moved very fast when the charge of the black hole is changed 
p8| . However the s = — 1 Kerr black hole modes move slowly even when the angular momentum of the black hole 
is changed. This can be understood, because the electromagnetic perturbation depends more on the change of the 
black hole electronic charge rather than the change of the angular momentum. Even though one might expect that 
the gravitational and the electromagnetic mode coalesce in the extreme limit as was established in the charged black 
hole case |^,^9|, obviously there is no such correlation between these two modes in Fig. 1. 

The behavior for the first few modes were rather simple, but the rapidly damped modes show the same strange 
behavior as in the Reissner-Nordstrom black hole. As presented in Fig. 2, the first strange feature appears for n = 3 
in the / = 1 electromagnetic modes and for n = 6 in the I — 2 gravitational modes. These modes go through a couple 
of small loops when the angular momentum increases. Similar loops are also found in the results for the higher mode 
number, n. In these figures it is clear that the loops are getting bigger and the number of cycles increases when the 
mode number is getting higher, as in the case of the Reissner-Nordstrom black hole. These strange features typically 
appear when the ratio of the damping rate to the oscillation frequency reaches a certain value, approximately between 
five and ten for any multipole index /, which is rather large compared with the charged black hole case. Since the 
highly damped modes are not well understood, it is difficult to explain why such a peculiar behavior occurs. In the 
previous paper p3] we suggested that the fall-off rate of the effective potential at the horizon might be the cause. This 
explanation suggests itself also in this case, because for the Kerr black hole the fall-off rate of the potential depends 
on the Kerr parameter in a similar way as the Reissner-Nordstrom black hole. 

It should be noted that we have tracked the trajectory of the I = 2 ninth overtone (n — 8) gravitational quasinormal 
mode in Fig. 2. This frequency is known to coalesce with the algebraically special mode at the Schwarzschild limit 
which is located on the imaginary axis. The algebraically special mode is defined in a completely different way from 
the quasinormal mode. It is defined as a mode that is purely ingoing at both boundaries. Thus , it se ems that there 
is no reason that these modes have any kind of relationship. We will discuss this matter in Sec. Ill C . 

The numerical error in the previous work 28 1 , possibly caused by the bad convergence property of the continued 
fraction, does not effect the Kerr black hole case. The convergence of the continued fraction is very fast for the Kerr 
black hole and thus can avoid numerical errors. When checking our code for the fundamental I = 2 gravitational 
mode, it was possible to obtain the mode accurately until a/M = 0.999995. For the Reissner-Nordstrom black hole 
the same mode could be calculated only up to Q/M = 0.9985. Thus, it is concluded that although the Teukolsky 
equations are more complicated than the Zerilli-Moncrief equations the quasinormal frequencies can be calculated 
very precisely through the continued fraction method even when the Kerr parameter is near the maximum value. The 
recurrence relations in the Reissner-Nordstrom black hole are four-term recurrence relations and numerical instability 
might be caused through the transformation of four-term relations into three-term relations. 



B. The cases of m / 

First we discuss the case of the s = — 2 gravitational perturbation with the least multipole index 1 — 2. ( The 
case with higher multipole index is similar. ) It is noted that m must be in a range of — / < to < I, thus we have 
four different modes for I = 2 with m ^= 0. Owing to the symmetry that we mentioned before, all the modes with 
the azimuthal multipole index m are obtained by reflecting the modes with the index —to about the imaginary axis. 
We therefore discuss only the modes in the right-half of the complex plane. The trajectories of the m = +1 modes 
are plotted in Fig. 3. We can find that the oscillation frequency becomes very fast when the Kerr parameter is 
increased, on the other hand, the damping rates become slower. The mode moves fast away from the corresponding 
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Schwarzschild mode. However the branch of to = — 1 moves slowly and stays near the corresponding Schwarzschild 
modes. It is surprising that the m = —1 mode turns sharply at some specific Kerr parameter approximately a = 0.2, 
which is not yet found for any non-rotating black hole. It is also noted that the mode of m = ±1 does not experience 
loops unlike a highly damped mode of the m = mode. This distinguishes the behavior of to 7^ modes from that 
of the to = modes. 

The behaviour of the m = ±2 modes is quite similar to that of the to — ±1 modes. Fig. 4 shows how these modes 
move as the Kerr parameter increases. As is easily seen, the to = +2 modes move very fast compared with any other 
modes and these modes converge to a specific undamped frequency at the Kerr limit. It is already known that the 
modes of the I — 2, m — +2 branch accumulate onto the critical frequency of super-radiance, oj c = 2. Detweiler has 
shown that when a = 1/2 an infinite number of modes cluster at the critical frequency. This suggests a marginally 



I, and this is quite interesting 
ose to the maximum the Kerr 



unstable feature of the extreme Kerr black hole discussed by several authors J3 1|— 133 
to investigate in detail. At present people believe that when the Kerr parameter is c 
parameter is decreased emitting the gravitational waves and losing rotational energy. Then the black hole settles 
down to a stable non-extreme Kerr black hole. You can also notice in Fig. 4 that not all the modes converge to the 
real critical frequency, but the sixth mode turns sharply at about a — 0.365 and then limits to a complex value around 
lu = (0.239, —0.242). Detweiler's result was that there are an infinite number of quasinormal modes converging to the 
critical frequency. Our numerical results suggest that most of modes converge into the critical frequency but some 
modes still stay away from the critical frequency. We so far have calculated the modes up to n — 20 but there was no 
other mode which does not converge into the critical frequency. 

It is also interesting that there seems to be some relation between the m = ±2 modes and the m = ±1 modes. 
As is shown in Fig. 5, the both to = +1 and to = +2 modes start from the Schwarzschild mode initially in the 
same direction and the both to = —1 and m = —2 modes moves in the opposite direction to the above two modes. 
Therefore the trajectories of these two modes are tangent to each other. On the other hand, the to = mode moves 
into a non-correlated direction. It seems that the to = ±2 modes move just twice as fast as the to = ±1 mode initially. 
(This can be also confirmed by examining the I = 3 case, where the m = 3 mode moves three times as fast as the 
to = 1 mode.) These feature can be understood, the m = ±2 modes are more distorted than the to = ±1 modes and 
hence the to = ±2 are more sensitive to the change of the Kerr parameter that corresponds how much the spacetime 
is twisted. 

Fig. 6 is a figure of the trajectories of the I = l,m = ±1 modes for the electromagnetic perturbations (s = —1). It is 
very similar to the case of the — s = to = 2 modes. These modes also converges into the undamped critical frequency 
u> c = 1. Owing to Detweiler's argument, for any s the I = m modes have a similar accumulation of an infinite number 
of the quasinormal modes onto the critical frequency u> = I. There was one mode which never converge into the 
critical frequency for the s = — 2 case, however for the s = — 1 case no such mode was found up to n = 8. 



C. Algebraically special modes 

Here we only consider the algebraically special mode of the s = —2 gravitational perturbations, because the only 
known implication that there is a correlation between the algebraically special mode and the quasinormal mode is 
that both modes for the Schwarzschild black hole coalesce for gravitation al p erturbations. 

The algebraically special perturbations were first considered by Wald |p4| , who was interested in how the pertur- 
bation can be the only independent component among several Weyl scalars. (Actually it is impossible if you consider 
only the real frequency.) As is well known, the perturbations have generally two independent components \&o and , 
which correspond to the ingoing and the outgoing wave, respectively. Thus the algebraically special mode is the mode 
that has only an ingoing or outgoing wave. Such kind of perturbations were first analytically solved by Chandrasekhar 
p0| . For the Schwarzschild black hole, algebraically special frequency can be easily solved by imposing the condition 
that Starobinsky's constant vanishes. (Because two modes are related to each other with a Starobinsky constant.) It 
is then found that there is such a mode on the imaginary axis at w = —2i/M. 

For the Kerr black hole it is rather difficult to obtain the algebraically special frequencies owing to the difficulty to 
get the separation constant Ai m without resorting to numerical computation. Nevertheless, if we use the continued 
fraction technique for the angular part of the Teukolsky equations, the separation constant is accurately obtained. 
The square of the Starobinsky constant is given by 

|Cf = A 2 (A + 2) 2 - 8uj 2 X{a 2 {5\ + 6) - 12a 2 ) + 1Uuj 2 {M 2 + a 2 + a 4 ), (3.1) 

where 



.,2 



am 



<\ = <r , (3-2) 
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and A, the separation constant in Chandrasekhar's notation, is related to ours with the relation 

A = Ai m — 2aujm + a 2 uj 2 . (3-3) 



Therefore we only need to solve \C\ 2 = together with Eq. (2.29). In contrast to quasinormal modes, algebraically 
special modes can be calculated even when the black hole is extreme, because we do not need to solve Eq. (2.3) which 
has two confluent singularities in the extremal limit. 

The solutions al ways appe ar as a pair, owing to a reflection symmetry about the real axis. It is also noticed that the 
symmetry of Eqs. (3.1) and ( 2.29 ) allows us to get the rn = —1, —2 modes by reflecting the m = +1, +2 modes about 
the imaginary axis as in the quasinormal mode. Therefore, it is enough to see the first quadrant. The solutions of I = 2 
branch which are obtained using the continued fraction scheme are plotted in Fig. 7 together with Chandrasekhar's 
result ]3C( ], and these values are listed in Table 1. The m = mode increases along the imaginary axis initially as the 
Kerr parameter increases and then goes away from the imaginary axis when a = 0.247222. In Ref. |3(i(| the value of 
the m = mode for a = 0.25 still stays on the imaginary axis, but Chandrasekhar mentioned that this mode does not 
satisfy a positivity of k. This contradiction was caused, because his values were obtained using a separation constant 
with large numerical error, and thus his value was wrong. Our values are all in agreement with this positivity check, 
since our scheme can avoid numerical errors. If you see Fig. 6, there is no m = -1 and m — —2 mode in the first 
quadrant, and the m — +1 and m = +2 modes go horizontally from the imaginary axis like the quasinormal mode. 
Despite that, there seems to be no correlation with any quasinormal mode. 

From the point of view of the scattering problem, a quasinormal mode is regarded as a pole of the reflection 
amplitude R(co) and the transmission amplitude T(lo), on the other hand, an algebraically special mode corresponds 
to a zero of R(ui). Therefore these two modes will never coalesce in general, which is consistent with the above 
numerical results. Superficially it might appear that these two modes coalesce at the Schwarzschild limit, but we 
believe that the quasinormal mode disappears at the limit owing to a cancellation with the corresponding mode in 
the left-half of the complex plane. 



IV. CONCLUDING REMARKS 



In this paper we have presented several hints to be useful for resolving a couple of questions concerning the 
quasinormal mode of the black hole. We have presented the results of detailed calculations of the quasinormal modes 
of the Kerr black hole. This supplements previous studies in several ways. First, we take the rapidly damped mode 
into consideration. Second, we also investigate precisely how these modes move as the Kerr parameter approaches the 
limit. Third, we also consider the relationship between the algebraically special mode and the quasinormal mode. 

As was the case for the Reissner-Nordstrom black hole p8| , the highly damped m = Kerr quasinormal modes 
generally spiral into the extreme values. Though not presented in this paper, we have found the similar behavior 
also in the case of the test scalar field on the Kerr background (s = 0). Therefore these spirals must be one of the 
fundamental features for rapidly damped modes. The behavior of the m ^ modes are, on the other hand, different 
from the behavior of the m = mode. The damping rates of < m < I modes decrease very fast without going 
through spirals. An infinite number of the m = I modes tend to accumulate onto the critical frequency, lu c = I, as 
a/M — > 1. We also found that there is one mode which does not converge to the frequency. Whether another mode 
like that exists is an open question. At present there is no evidence up to n = 20 for I = 2 gravitational perturbations. 
The — I < m < modes have rather simple behavior in contrast to the positive m cases. These modes move slowly 
around the corresponding Schwarzschild modes. A interesting fact is that some of the modes go through a quick turn. 
These surprising features for ra/0 were not yet found for the spherical black holes and hence should be considered 
in detail. Since the m ^ modes do not exist for non-rotating black holes (all the m ^ modes degenerate into 
the m = case owing to the spherical symmetry), these modes are actually specific for the rotating black holes. We 
also obtained the precise values of the algebraically special modes using a continued fraction method. There was no 
correlation between the algebraically special mode and the quasinormal mode except that these modes superficially 
coalesce in the Schwarzschild limit. Currently we do not know why the coalescence occurs. However, owing to the 
existence of the algebraically special mode at lo = — 4i the quasinormal mode can not be there. Two quasinormal 
modes coalesce and then probably cancel at the limit. This cancellation prevents the existence of the purely damped 
quasinormal mode. 
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FIG. 1. Solid lines and dashed lines are trajectories of the m=0 quasinormal frequencies of the gravitational and elec- 
tromagnetic wave, respectively. Each right endpoint corresponds to the Schwarzschild quasinormal frequency, and each right 
endpoint corresponds to the frequency in the limit of nearly maximal angular momentum a/M=0.99. Ticks marked on each 
line are the frequencies for a=0.1, 0.2, 0.3 and 0.4 from the left to right. 
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FIG. 2. The first three specific quasinormal modes for m=0 are plotted. The marks correspond to the frequencies of 
a=0, 0.05, 0.10, 0.45. These modes are generally spiral into the values in the extreme limit. The qualitative features for 
the gravitational and electromagnetic frequencies are basically identical. These figures are much finer than the figures in the 
previous work for the charged black hole because the continued fraction scheme works fine with less numerical errors. 
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FIG. 3. The first nine m=-l and m=+l gravitational quasinormal modes for 1=2 are plotted. The diamonds are the 
Schwarzschild quasinormal frequencies and the right part starting from the Schwarzschild frequency is the branch of the m=+l 
modes, which is depicted by the dashed lines. The left part is the branch of the m=-l, which is depicted by the solid lines. 
We marked the values of the a=0.05, 0.1, 0.15, 0.45 Kerr frequencies for the m=+l branch and a=0.1, 0.2, 0.3, 0.4 for the 
m=+l branch. It is noticed that the ninth mode comes out of the imaginary axis that corresponds to the algebraically special 
mode of the Schwarzschild black hole. 
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FIG. 4. The first nine 1=2 and m=-2,+2 gravitational quasinormal modes are plotted. All the modes of the m=+2 branch 
except the sixth overtone (n=5) converge into the critical frequency on the real axis. This accumulation of modes is in agreement 
with Detweiler's result. But our numerical results suggest that some of the modes do not cluster into the critical frequency. 
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FIG. 6. The first nine m=-l and m=+l quasinormal modes for the electromagnetic waves in the Kerr geometry are shown. 
This is very similar to the case of l=m=2 gravitational quasinormal modes, however all the mode up to n=8 converge into the 
critical frequency in contrast to the gravitational frequencies. 
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FIG. 7. The algebraically special modes of m=0, +1 and +2 for 1=2. There also exist another branches of modes symmetric 
about the real axis. The mode of m=-l and m=-2 are in the left half plane, which are obtained by reflecting the m=+l and 
m=+2 modes about the imaginary axis, respectively. The m=0 algebraically special mode first goes up vertically and then 
produces two branches, one in the right half of the complex plane and the other in the left half of the complex plane, which is 
not shown in this figure. The corresponding Kerr parameter is a=0. 247222 and the pure imaginary frequency is about -6.681L 
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(1.7626, -2.4766) 




0.35 


(2.84578, -4.36141) 


(1.71091,-2.18596) 


(1.41732,-1.50403) 






(1.53041,-1.9104) 




0.40 


(2.93170, -3.71616) 


(1.69372,-1.97997) 


(1.36428,-1.36084) 


0.45 


(2.91142, -3.23046) 


(1.66300,-1.80540) 


(1.31270,-1.24256) 


0.50 


(2.84632, -2.85259) 


(1.62505,-1.65663) 


(1.26369,-1.14328) 



TABLE I. The algebraically special modes calculated through the continued fraction algorithm are compared with those 
obtained in Chandrasekhar's paper. The upper values are obtained through the continued fraction method and the lower 
values are from his paper. It seems his values are not accurate for large Kerr parameters, owing to the numerical errors caused 
when getting the separation constants. 
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